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Abstract 

A method allowing analysis of the overpopulation of phase-space in heavy ion collisions in a 
model independent way is proposed within the hydrodynamic approach. It makes it possible to 
extract a chemical potential of thermal pions at freeze out irrespective of the form of freeze-out 
(isothermal) hypersurface in Minkowski space and transverse flows on it. The contributions 
of resonance (with masses up to 2 GeV) decays to spectra, interferometry volumes and phase- 
space densities are calculated and discussed in detail. The estimates of average phase-space 
densities and chemical potentials of thermal pions are obtained for SPS and RHIC energies. 
They demonstrate that multibosonic phenomena at those energies might be considered as a 
correction factor rather than as a significant physical effect. The analysis of the evolution of the 
pion average phase-space density in chemically frozen hadron systems shows that it is almost 
constant or slightly increases with time while the particle density and phase-space density at 
each space point drops down rapidly during the system's expansion. We found that, unlike the 
particle density, the average phase-space density has no direct link to the freeze-out criterion 
and final thermodynamic parameters, being connected rather to the initial phase-space density 
of hadronic matter formed in relativistic nucleus- nucleus collisions. 
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1 Introduction 

The interferometry method to measure particle phase space densities averaged over the volume of 
a system that emits the particles was first proposed by Bertsch If systems formed in ultra- 
relativistic A+A collisions are thermal, then it is possible to compare the experimental value of 
the average phase-space density (APSD) with one associated with the Bose-Einstein 2^ or Fermi- 
Dirac pj thermal distributions, /be,fd(p), and extract in this way a particle chemical potential 
(or fugacity) at the freeze-out temperature. Such a chemical potential could have a positive value 
even for pions since a hadron gas evolution is chemically frozen rather than chemically equilibrated 
@J El HI Ej • If it turns out that the particle fugacity is fairly large, then the boson density could 
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be around the critical one and significant multibosonic I1U1 1111 112j . and coherent ^3] effects 
might be present at low particle momenta. 

There are, however, two serious problems which have to be solved to apply the above method. 
It is obvious that the APSD (f(p)) can be associated directly with a thermally equilibrated distri- 
bution for static homogeneous fireballs, /be(p) for bosons. For expanding inhomogeneous systems, 
characterized by a set of (local) equilibrium distributions which are Doppler shifted in the vicinities 
of each space-time point according to the collective velocity field, the phase-space density averaged 
over some hypersurface situated at or "after" the freeze-out could lose, in general, the direct con- 
nection to the Bose-Einstein or Fermi-Dirac phase-space distributions. As it was demonstrated in 
Ref. j!4j . even at a uniform temperature and chemical potential at freeze-out, the intensity and 
profile of longitudinal and transverse flows influence very essentially the value of {f(p)}, making, 
thus, the extraction of fugacity by comparison with fBE(p) ambiguous and model dependent. One 
of the aims of this paper is to study and substantiate the model independent possibility for an 
analysis of overpopulation of the phase-space in 3D expanding fireballs at freeze-out. 

The another problem is resonance contributions to particle spectra and interferometry radii 
15 . If chemical freeze-out takes place, then about 2/3 of pions comes from decays of resonances in 
high energy heavy ion collisions El • Only a small part of them are long-lived resonances which 
give the effect of a halo ^0] - a suppression of correlation functions of identical pions (see also ) • 
The others, short-lived, enhance significantly the pion single-particle spectrum, as compare to the 
thermal pion one, and also modify the interferometry radii. Our analysis of pion production by 
some concrete resonance species shows that even for thermal distribution of those resonance species 
at freeze-out the spectrum of the produced pions is far from the thermal one. Also, an effective 
region of such a pion emission at small momenta is not associated with the interferometry volume 
at the corresponding small mean momenta of pion pairs from the resonance decays. Peculiarities 
of pion production by resonances are important to understand the influence of resonance decays on 
the average phase-space densities which are currently measured in experiments at CERN SPS and 
BNL RHIC. We analyse quantitatively the total effect on spectra and interferometry volume of pion 
production by huge number of resonances with masses up to 2 GeV and dependence of the effect 
on the intensity of flows. As a result, we propose how to cope with the resonance contributions to 
study the freeze-out properties of the expanding systems. 

One of the most serious challenges to the theoretical understanding of the processes of A+A 
collisions is multiplicity (and energy) dependence of the APSD (/) in central events. According to 
the Bertsch result 1 (/) ~ Jv^f? ■ If at different energies of A+A collisions the pion interferometry 

volume Vi n t of the created systems is proportional to the pion multiplicity dN/dy, and T e ff, the 
effective temperature of the transverse spectra of pions, changes only slightly, then approximate 
constancy of the APSD could be considered a the natural property of the system at freeze-out in 
high energy A+A collisions. It was a very likely hypothesis at relatively small energies |2j, since 
it had been supported by the observed proportionally between dN/dy and Vi n t JZj. However, this 
proportionally law starts to be violated already at SPS energies and is completely destroyed in 
experiments at RHIC |18j . In the latter case the interferometry radii and volume turn out to be 
essentially smaller than was expected at the corresponding RHIC multiplicities (the HBT puzzle) 
19 . As a consequence, the APSD at RHIC energies is found to be the considerably higher than 
at the formerly achieved energies |2()j . Does this mean that freeze-out conditions of the system in 
heavy ion collisions at RHIC are essentially different from the ones at SPS? Does this signify that 
the freeze-out particle densities at RHIC come close to the critical value? The latter could lead to 
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the multibosonic effects at freeze-out which, in turn, can reduce the observed HBT radii [HI HUj. 
Are thus multibosonic effects relevant to the HBT puzzle? 

To understand the above problems we study the APSD and particle density temporal evolution 
in a few hydrodynamic models. We observe, in particular, that the phase-space density averaged 
over the freeze-out volume, (f(p)), drops rather slowly during the 3D expansion of the system, - 
much slower than the phase-space density f(t, r, p) at any fixed r and than the particle density. 
This difference in the behavior of the two groups of values becomes bigger when the intensity of 
the transverse flow grows. The detailed analysis of this effect as well as the model independent 
estimates of chemical potential for the thermal pions at freeze-out makes it possible to shed light 
on the reason of the relatively high APSD at RHIC and on the HBT puzzle. 



2 Evolution of the average phase-space density. 

2.1 Temporal behaviour of APSD (f(p)) in analytically solved hydrodynamic 
models 

Let us demonstrate the important features of the APSD evolution versus particle density one by 
one using hydrodynamic models of A+A collisions. To see the main idea we start from the exact 
solutions of perfect non-relativistic hydrodynamics. If the initial conditions correspond to the 
thermal Boltzmann distribution with uniform temperature To, zero initial flows, u(t = 0, r) = 0, and 
spherically symmetric Gaussian (with radius i?o) profile of particle density, then the hydrodynamic 
evolution of that system is described by the local equilibrium phase-space distribution which is an 
exact solution of Boltzmann equation for any cross-section of particle interactions [U and has the 
form: 



R 3 \(2TT) 2 mT J ^ 2T(t) T(ty 
where N is total number of particles, v = p/m are their velocities, and T{t) = 1+To Jy mR $ > f^y = 

r 2 mRl mfigx , n _ tT(t) 

One can check directly that the momentum spectrum and the interferometry radii, evaluated at 
any time t (at any isotherm) are identical to those calculated at the initial time t = 0. The easiest 
way to see it is to rewrite expression (0) in the form 
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where the effective temperature of the spectra is a constant in time and equal to the initial one, 
T e ff = Tq. Thus, the maximal phase space density and Bertsch's average phase-space density 

is the same at each p as it was at the initial time! 

The space density of the particles falls down with time as 1/i 3 at large t. So, the particles in 
the system will really stop interacting at some later stage and the system decouples, preserving its 
maximal PSD and average PSD values as they were at the moment when the fireball was created. 
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It is noteworthy that the local phase-space density f Leq '(t, r, p) in the system drops with time at 
any fixed point r very rapidly, as exp[— (c 2 t 2 + C2t)]. 

The perfect constancy of the APSD is caused by the peculiarities of this specific model: in 
particular, its spherical symmetry, as was discussed in detail in Ref. |22| . At the same time, as we 
will demonstrate, the important feature of this simple model, namely, the principal distinction in 
temporal behaviours of different types of densities, is preserved in realistic models of A+A collisions: 
the APSD decreases with time essentially slower than the particle density. 

To make the next step towards clarifying the effect, let us consider one more analytical il- 
lustration of it based on the model of a 3D expanding fireball with longitudinally Bjorken-like 
flow v z = r z /t . Here is the exact non-relativistic solution of the ideal hydro equations for an 
initially Gaussian transverse density profile which results in the following evolution of the locally 
equilibriated phase-space density [SJ: 



f Leq -(t,r,p) 
n(t, r) 
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where the initial particle density in the central part of system, r x = r y = 0, is no = N / (2tt) 3 ^ 2 R^tQ 
and 

( figto N 2/3 - - T 



At each moment t the longitudinal spectrum in this model is m^- = 2irR I ) tonQ = const, similar to 

dN 
dy 



^ = const in the relativistic boost-invariant model, and the momentum spectrum is 



d 3 N n Rlt . V T x (R \ 
- 1 3- = ^= — --exp -— — — — , 6) 
dtp m z T e ff(t) 2mT e ff(t) 

where the effective temperature, T e ff, is 

T eff = m R +T. (7) 

The T e ff(t) does not change much since the decrease of temperature T when the system expands 

is accompanied by the increase of transverse flows described by R in Eq. (JHJ). The typical behavior 
of the temperature and effective temperature in this model can be seen in Fig. 1. 
The interferometry radii, out-, side- and long- take the forms 

R o = R s = — R > R L = t — ( 8 ) 

l ef f m 

and are similar to the ones in the relativistic case at the central rapidity point [231 124j , if m — > 



mr = ym 2 +Px- The interferometry volume, taking into account Q, is 

R 2 MT a f/ 2 



V int = RoRsRL = " UV (9) 
T effVrn 
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and does not change significantly since it depends on T e ff only. So, it is not surprising after all, 
that the APSD is, at least approximately, preserved with time 



f 1 \ 3/2 o 2 1 

</(*,p)> = n l-T—r) ex P(-?^7T)- ( 10 ) 
\4irmToJ 2mT eff (t) 

The average momentum phase-space density, (f(t,p)}, as well as the interferometry volume Vi n t, 
will be preserved in time if and when, at some stage of the evolution, T e ff(t) is constant in time. 

Practically, the effective temperature, and thus the APSD (f(t, p)), could be slightly reduced or 
even increased with time during the evolution, depending on initial conditions, the evolution stage 
and particle mass. Behavior of the effective temperature is correlated unambiguously with the be- 
havior of particle transverse energy per unit of rapidity (or average particle transverse momentum) . 
In the idealized case of the ID longitudinal expansion of a one-component gas, the particle trans- 
verse energy per unit of rapidity can only decrease since the effective temperature of the transverse 
spectrum in that case is just the system's temperature that falls down permanently. However, in 
a realistic case when transverse expansion is switched on, the transformation of a system's heat 
energy into the transverse and longitudinal kinetic energies of particles becomes much more com- 
plicated. In the 3D expansion the diminution of the total energy of the part of the hydrodynamic 
tube confined within a small rapidity interval is accompanied by a faster decrease in the particle 
longitudinal kinetic energy, comparing to the ID case. It is caused by faster temperature decrease, 
while longitudinal collective velocities are the same in that piece of fluid. As for the transverse 
energy, it accounts for not only a decrease of the heat energy, associated with temperature, but 
also an increase with time of transverse flows. Therefore, the particle transverse energy (in the 
relativistic case it is better to consider average transverse momentum of particles) essentially does 
not change during the evolution. Such a process of transformation of heat into transverse flow 
is accompanied by an appearance of positive correlations between the absolute value of particle 
transverse momentum and its distance from the centre of the system. The degree of that corre- 
lation grows with time compensating for the contribution of increasing geometrical radii to the 
interferometry volume Vi n t(pT) Since the transverse momentum spectra do not change much, 
the APSD (f(t,p)) ~ dydtpT l^intipr) becomes about constant at some stage of evolution unlike 
the local phase-space density f(t, r, p) and particle density. The latter is just a fast monotonically 
decreasing function (see Q)): 

This is one of our main points - the system can reach small densities and decay with rather high 
APSD reflecting the high average phase-space density at the moment when the system was formed! 

Now we pass from general physical ideas to the concrete results within the realistic model of 
Pb+Pb collisions at CERN energies 158 AGeV [Sj. The model describes hydrodynamically the 
evolution of chemically frozen hadron-resonance gas from the hadronization stage to the thermal, 
or kinetic, freeze-out. As was found recently [3], if chemical freeze-out is incorporated into hydrody- 
namics, then the final spectra and fireball lifetimes are insensitive to the temperature at which the 
switch from hydrodynamics to cascade RQMD is made. This implies that the local conservation 
law describing the transformation of heat energy into collective flows for a chemically frozen hydro- 
dynamic evolution results in spectra which correspond approximately to the microscopic cascade 
calculations. Based on an analysis of the particle number ratios, pion rapidity density at y ~ 0, 
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slopes of transverse spectra and values of interferometry radii for pions, kaons and protons, the 
reconstruction of the hadronization stage was done in Ref. JB] by solving backwards in time the 
hydrodynamic equations for the chemically frozen composition of mesons and baryons with masses 
up to 2 GeV. The "initial conditions" for such a solution were defined from an analysis of spectra 
and correlations at the thermal freeze-out stage. Supposing that the hadronization of Bjorken-like 
expanding quark-gluon plasma (QGP) happens at a uniform chemical freeze-out temperature T c ^ 
with Gaussian-like particle density transverse profile of some width R and self-similar transverse 

flow ut = r R /R , one can restore from the hydrodynamic solutions the hadronization proper time 
T c h, transverse radius R(r c h) , transverse velocity VR{r c h)=R (r c h) at the Gaussian "boundary" of 
the system, and then the energy and particle number densities in the center of the system, e c h and 
n c h- The temperature T c h and baryonic chemical potential HB(T c h) are defined in Ref. (Hj from the 
analysis of particle number ratios. 

TABLE 1. The characteristics of the chemical freeze-out in Pb+Pb SPS collisions at 158 AGeV 

taken from Ref. jHj. 



T ch GeV 


fJ-B,ch GeV 


T c h fm 


Rch fm 


VR,ch 


_ GeV 
t-ch J^J 


Tlch 7^3" 


0.164 


0.224 


7.24 


4.63 


0.304 


0.420 


0.421 



The evolution of the APSD {f(r,p)) for the frozen number of direct pions in hadron-resonance 
gas with the initial conditions at chemical freeze-out as in Table 1 is demonstrated in Fig. 2. 

As one can see, the APSD is changed less than 10 per cent during the evolution from the 
hadronization stage to the kinetic freeze-out and its values at small transverse momenta are even 
higher at the final than at the initial moment. At the same time, as is shown in Fig. 3, the 
densities of direct pions drop rather quickly in the central part of the system where soft momentum 
particles are produced. One can see the principal difference in the evolutions of the space- and 
average phase-space densities from the comparison of correspondent curves at r = 15 fm/c and 
r = 7.24 fm/c. While the APSD changes not more than 35 per cent and even increases at small 
Pt, the particle density drops a factor of 6 ! We can conclude, therefore, that the basic physical 
phenomenon of approximate "conservation" of the pion APSD during the system evolution that 
we demonstrated analytically in the non-relativistic models of the 3D hydrodynamic evolution also 
takes place in the general case of relativistic chemically frozen 3D expansion of a many-component 
hadron resonance gas. 



2.2 Totally averaged phase-space density (/) as integral of motion 

Let us devote the end of this Section to the new, important results concerning the evolution of 
the phase-space density {f(t)) that is totally averaged over momentum and space with the integral 
measure d 3 ra?p : 

- I f 2 (t,r,p)d 3 rd 3 p 
U{ )) - Jf(t,r,p)d 3 rd 3 p' [U) 

It is obvious that (f(t)) ^ f d 3 p {f(t, p)) . Supposing the Boltzmann non-relativistic approximation 
for the distribution function holds: 

/( t ,r,p) = (2 J )-3exp( ro < V m y' )^r) (13) 
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where u(t, r) is the collective velocity field, £(t,r) = exp(/x/T) is the fugacity, and \i is the (non- 
relativistic) chemical potential, one can find for the finite expanding system that 4f (f(t)) = 0. 

Indeed, from the non-relativistic hydrodynamic equations for the one-component Boltzmann 
gas (see, e.g., (23]) it follows the fugacity £ is preserved along the current lines: (dt + ify,d ri )£ = 0. 
Then one should appeal again to the hydro equations to deduce that 

(d t + d nUi )T 3 / 2 ^ = 0, (14) 

and get, finally, using the property of finiteness of the system, that (k = 1,2, ...) 

/ d3r ¥t (/ •f*fr r 'P) rf3 f) = 

(15) 

So the totally averaged phase-space density is conserved during the system expansion! Such models 
which describe initially finite systems, f(to, r — > oo) — ► 0, belong to the class of Landau-like models 
26 . Of course, the property {fit)) = const is satisfied in the concrete model ^ of spherically 
expanding fireball. 

There is, however, another, Bjorken-like class of hydrodynamic models, which are boost invari- 
ant and, so, formally infinite in the longitudinal direction. For this situation it is better to define 
the phase-space densities averaged over all phase-space except particle longitudinal momentum p z 
, or rapidity y, in order to associate the APSD with the small rapidity interval near mid-rapidity, 
vl = Pz/m ~ 0, or y ~ 0. We will mark the corresponding APSD as (f(t,y)) . Since the boost- 
invariance dictates that distribution function depends just on the difference vj_, — (ul = r z /t) 
for the non-relativistic case, or y — n (rj is hydrodynamic rapidity ) for the relativistic one, the 
integration over fluid longitudinal coordinates at fixed particle velocities vl ~ 0, or y 0, in (fT2l 
can be substituted with the integration over particle rapidity at fixed hydrodynamic velocities ul 
~ 0, or n ~ 0. Therefore, (f(t,y)) corresponds to a totally averaged phase-space density in a small 
central slice An of the hydrodynamic tube that expands conserving entropy S and particle number 
(dN/dn)An. Using Eq. (|14[) one can find, similar to ()15|) . that 

td 2 r J f k (t,r,p)d 3 p^J=0, (16) 

and, thus, the APSD (f(t,y)) in this case is also preserved in time. 

It is useful to explain and generalize to the relativistic case the above results without resorting to 
hydrodynamic equations. Let us consider the time evolution of the set of i— fluid elements starting 
from some initial time t$. They contain particle numbers AAj, entropy AS{, energy in their rest 
frames (marked by asterisk) AE*^ and occupy initially the small volumes AV*, \ near points Xj( ) 
with temperatures Tj(o) and chemical potentials /%(o). At some later time tf the fluid elements are 
characterized by the corresponding values: (0) — > (/) while preserving their entropy ASi (due to 
isentropic expansion) and particle number ANi (because of chemically frozen evolution). Then, 
according to the thermodynamic identity (p is pressure) 

^ = PwAVfo + Agp) _ m ^ = P^AVfo + AE^ _ 

T m T i{0) T m T i{f) T i[f) T i[f) 



d_ 

It 



d 3 r / f k (t,r,p)d 3 p 
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In the ideal Boltzmann gas p^* = NT and at m/T 3> 1 E « mAf + liVT. In the ultrarelativistic 
case m/T< 1 the particle mass can be neglected and i£ ~ 3NT. As the result of Eq. I|17|) j m 
both cases the " non-relativistic fugacities" will be conserved, m = m , along the world 
lines of each fluid element. What happens to the APSD during the system evolution? In the 
case of relativistic expansion the averaging in (j!2j) should be done in relativistic covariant form 
cftrcPp^-da^p^^-^ , f(t) — ► f(o~), f(t,vi) — > f{o~,y) where the hypersurface a generalizes the 
hypersurface of constant time t. The 4-vector formed by the world lines of the i— fluid element 
crossing a and moving with 4- velocity is f^y* da^ = , where n^(x) is normal to o~(x). 

Then (k = 1,2,...) 

/ da" / pj^- e W [-k(p^ - m)/T]i k (18) 
Jav* J Po 

~ Tt /2 £ k AV? if m/T > 1 and ~ i; 3 £ fc Ay.* if m/T <C 1. Since the "non-relativistic fugacity" £ = 
expf^^^] and corresponding values (mTi) 3 / 2 AV* ( m/T 3> 1) or AT^* (m/T <C 1) are constants 
along the current line, then 



d <o) / Pm— / = / d*J / P„— f (19) 

and, finally, the values of phase-space densities (/(<r)) and (f(a,y)) averaged over do^-^ are 
the same on different hypersurfaces. In other words, the totally averaged PSD of the relativistically 
expanding Boltzmann gas is approximately preserved during the system's evolution if the mass of 
the particles is much higher or smaller than the temperature of the system. 

The non-relativistic Bjorken-like model (@J) gives an example of such a conservation of the 
APSD, resulting in m (f(t, vl)) = 8 v ^F(mT ) 3 / 2 = cons ^- The Pi° n APSD (f(t,y)) in the relativistic 
chemically frozen 3D Bjorken model 8, with the initial conditions as in Table 1 is also roughly 
constant: it changes from the chemical to kinetic freeze-out stages 5 — 7 per cent only, but it 
increases with time even without feeddown from resonance decays! This means that the non- 
relativistic fugacity of pions, expf ^ ^ ri7r ], grows with time. This paradoxical behavior of the APSD 
is, however, typical for a mixture of heavy and light hadron gases. Indeed, let us imagine that we 
have a thermal mixture of two gases: massless (e.g., neutrino) with concentration k\ and massive 
(e.g., neutrons) with concentration K2- K\ + K2 = 1. Each of them will preserve its non-relativistic 
fugacity if they expand separately. For evolution of the mixture we have, generalizing (|17|). that 
Kifii/T + ^2(^2 — m)/T = const along the current line at isentropic expansion. One could suppose 
that the expansion is conditioned mainly by the heavy component: (mTj) 3 / 2 AV^* ~ const during 
the evolution of the i-th fluid element. If then both fugacities are preserved then the entropy of the 

mixture will decrease: the light component needs the volume = A V(q) 7^- to preserve entropy 



T L AV? r 



at the same fugacity, but it gets only AV^) = A V(q) 7$p- At/ (0> < AV7~, while the heavy component 

preserves its entropy. Thus, the fugacity of the light component should grow and, correspondingly, 
the fugacity of the heavy component should decrease. The expansion of the mixture then follows 
a law which is intermediate between V ~ I/T 3 / 2 and V ~ 1/T 3 . Such behavior was considered 
in detail in Ref. (Hj. The mean mass of the chemically frozen hadron-resonance Boltzmann gas is 
significantly higher than the temperature of hadronization: m = ^ Kjirij = 0.662 GeV, where Kj 
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is the concentrations of j— hadronic species. The aforementioned result on the increase of the pion 
APSD (f(t,y)) during the chemically frozen hydrodynamic evolution is explained, therefore, by 
the fact that the pion gas is the lightest component in the mixture and so the pion non-relativistic 
fugacity should grow regardless of inelastic rescatterings and resonance decays. 

3 Resonance influence on spectra and interferometry volume in 
A+A collisions 

3.1 Emission and decays of resonances in quasi-classical approximation 

The above analysis is based on thermal local equilibrium distributions which in hydrodynamic 
models describe so called direct particles and resonances. The pion spectra from resonance decays 
are not thermal and, even under an assumption of sharp kinetic freeze-out, are formed in some ef- 
fective post freeze-out 4-volume that depends also on resonance widths. Unfortunately, one cannot 
ignore a priori the contributions to the pion spectra and interferometry radii by the short-lived 
resonances as they give more than half of total pion yields in ultrarelativistic A+A collisions. Such 
a big contribution arises since the chemically equilibriated composition of the hadron resonance 
gas, corresponding to relatively high hadronization temperature, is nearly frozen during the subse- 
quent evolution to low temperatures. Of course, the resonance decays could take place during the 
evolution of the dense system and lead to some increase of number of (quasi) stable particles even 
before the kinetic freeze-out. On the other hand, there are back reactions that produce the reso- 
nances in the dense system. These reactions do not lead to chemical equilibrium, but just reduce 
the effect of resonance decays into chemically non-equilibriated media. According to the results of 
Ref. [S], nearly 2/3 of all pions arise from the resonance decays after the thermal freeze-out if one 
totally ignores the reduction of resonance numbers during the hydro evolution. Taking into account 
a rather short duration of the hydrodynamic stage for the hadron gas, about 1.7 fm/c, which is 
comparable with the average life-time of the resonances, we can very roughly estimate that at least 
half of the pions arise from resonance decays during the post hydrodymanic stage. 

The only attempt to analyse the contribution of large numbers of short-lived resonances to the 
spectra and interferometry volume was done in Refs. [2ZJ within some hydrodynamic parametriza- 
tion of the kinetic freeze-out. The aim of those papers was to fit the spectra and interferometry data 
from CERN SPS and BNL RHIC supposing that kinetic freeze-out happens just after hadronization. 
However, the analysis of relative resonance contribution to the interferometry radii has not been 
done probably since a non-physical approximation corresponding to infinite widths (zero lifetimes) 
of all resonances was used. But such estimates are necessary for an analysis of the phase-space and 
particle density properties at kinetic freeze-out. Pursuing this aim we analyse the resonance contri- 
butions to the spectra and interferometry radii based on the hydrodynamic model of Ref. [Sj and 
take the widths T of mesonic and baryonic resonances with masses up to 2 Gev from Particle Data 
Group [2S]. Because of technical reasons we ignore, however, the contributions from many-cascade 
decays. We hope these do not change the general picture and conclusions. As we show below, the 
single particle transverse pion spectra have the same main features as in the model accounting for 
the many-cascade decays of heavy resonances |27j . 

To evaluate the contribution from resonance decays to the pion spectra and interferometry 
radii we use the quasi-classical approximation |291 13(J| for the emission function gi(x,p) of the i-th 
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resonance species: 



9i(x,p) = J du?Lpi(u?) j ^ 1 2m i 5((p i - pf - uj 2 ) ■ 

(20) 

oo 

J dr*r i exp(-r i r*) j da^)^ (x* - [xf + ^-\)f\ ^Pi) 

cr 

The emission function gt(x,p) describes the production of pions with momentum p from the i- 
th resonance species, the latter being distributed according to the local equilibrium phase-space 
density f\' eq '{xi,pi). The quasi-classical picture supposes that the i-th resonance is created from 
the hydrodynamic tube at a space-time point xf of the kinetic freeze-out hypersurface a and decays 
into tt + X after some proper time with the mean value 1/IV The probability to produce a group of 
particles X with an invariant mass u which accompany the pion at the decay of the i-th resonance 
is described by the covariant amplitude of the resonance decay ifi(uj 2 ) that is supposedly known. 
The pion spectra and correlations are defined by the phase-space densities of thermal ("direct") 
pions and resonances at a. According to Ref. [S], the local Boltzmann distribution j \' eq '(xi,pi) has 
a transverse Gaussian- like density profile with width R(rth), and is attributed to the kinetic freeze- 
out hypersurface a: Tth = const where the system is characterized by the constant temperature 
Tth, boost-invariant longitudinal flow and transverse flow i/t = vrt /R(jth)- The numerical values 
of correspondent parameters for Pb+Pb 158 AGeV collisions were extracted in Ref. [5] from an 
analysis of the spectra and interferometry radii at fairly large px to minimize the influence of 
resonance decays. They are T t h = 0.135 GeV, R(r t h) = 5.3 fm, vr = 0.4. We take proper freeze-out 
time Tth = 8.2 fm/c for the "direct" pions B., and use this freeze-out time also for resonances which 
produce pions. Note, however, that the average freeze-out time for different hadronic species was 
found to be 8.9 fm/c jH] which does not differ much from what we use. We do not intend in this paper 
to perform a detailed fitting of the spectra and radii; our aim is to analyse the dependence on flow 
intensity of the relative contributions of resonances to the spectra and interferometry radii within 
the appropriate model |Hj and, thus, to estimate the influence of resonances on the phase-space 
densities. So, we just take the parameters which were found in [S], and use them for highest SPS 
energies. Also, we estimate roughly the correspondent contribution to the APSD at RHIC energies 
by using the same model, simply enhancing the intensity of transverse flows: vr = 0.4 — » vr = 0.6. 



3.2 The results on spectra and interferometry and their interpretations 

The analysis of the single pion spectra from decays of the fixed resonance species shows it is es- 
sentially non-thermal, non-exponential, at least in the transverse momentum region below pq. One 
can see it at the top of Fig. 4 for the fixed heavy resonance species /2(1270) with po = 0.622 GeV. 
Nevertheless, the total contribution (thermal + decays) from many such resonances with different 
masses and different po results in a slope of the transverse spectrum which, as one can see from 
Fig. 4, is surprisingly similar to the spectrum slope for the thermal ("direct") pions in the wide 
momentum region at different intensities of transverse flows. The same peculiarities of the reso- 
nance contributions to the transverse pion spectra were found also in Ref. within the Landau 
hydrodynamic model with realistic freeze-out conditions. In contrast to these observations, the 
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results of Refs. |27| demonstrate a noticeable decrease of the effective temperature, or inverse of 
the slope of total spectra, as compared to thermal pions. The main reason for such a discrepancy 
is the (unrealistically) high decoupling temperature T that is taken to be equal to the hadroniza- 
tion, or chemical freeze-out one in Refs. |27| . It is qualitatively clear that the ratio of effective 
temperatures of total to thermal pion spectra depends on the relation between the mean po and 
T. We demonstrate this effect in Fig. 5 where one can see that the slopes of the total pion spectra 
are about the same as they are for thermal pions in the realistic interval of freeze-out temperatures 
100 — 140 MeV. Above this region the pion effective temperatures for the total transverse spectra 
are significantly smaller than for just thermal spectra. The result corresponding to the decoupling 
temperature equal to the chemical freeze-out one, 163 MeV [H], is close to what was found in 
Refs. [22j for similarly high temperatures T ~ 165 MeV. Clearly, at freeze-out temperatures below 
100 — 140 MeV the transverse effective temperature of total pion production has to be noticeably 
larger than in the case of thermal pions only. The correspondent tendency one can perceive from 
the plots in Fig. 5. 

Let us move to a discussion on the influence of resonance decays on Bose-Einstein correlations. 
Fig. 6 demonstrates that resonance decays extend the interferometry volume, especially at small 
transverse momenta. The ratio of the complete volume to the direct one grows when the intensity 
of transverse flow increases. It could be the double-effect of an increase in the resonance velocities 
in the cm. of the system and correspondent Lorentz dilation for the resonance life-times; both 
effects increase the mean free path for outgoing resonances. Another impressive effect of the flow 
and resonance decays on the interferometry radii is demonstrated in Fig. 7. One can see that the 
"standard" hydrodynamic behavior of the ratio Rout/Rside which is typically larger than unity 1 
becomes essentially smaller, R ut/Rside < 1> outside the region of small px if the resonance decays 
are taken into account. The above effects are created by the rather complicated interplay between 
the kinematics of resonance decays and flows. A description and understanding of the correspondent 
effects are crucial for the correct estimates of the phase-phase densities. 

Let us consider the effect of the relativistic longitudinal flow. If there is no flow and the 
resonance momentum spectra are, for instance, pure thermal, then all the points of the whole 
source volume are accessible with equal probability for resonance emission. Since the resonances 
have finite average life-times, they can decay inside as well as outside of the source volume, thereby, 
at any momenta the pion interferometry radii will be larger than the size of the source that emits 
the resonances. If the source expands relativistically in the longitudinal direction, then the spectra 
of heavy resonances are defined mainly by the flow, not the temperature. Let us suppose that 
pairs of identical pions with average rapidity y ~ and half-sum of pion transverse momenta pt 
are produced by two resonances 1* and 2* having, for simplicity, only two- particle modes of decay 
(pion plus, e.g., nucleon) with the same absolute value of pion momenta po in the rest frame of each 
resonance. If one studies the /ong-projection of the correlation function, the transverse momenta of 
both pions is supposed to equal pt • Since only pions with close momenta interfere, the resonances 



should have rapidities y* ~ ±y% = ± sinh y J± \m 2 ^° g uaran t ee 2/ — 0, if the resonance masses are 

much bigger than the temperature and there are no transverse flows. Thus, the heavy resonances 
with y\ = ±2/2 are emitted just from the two fluid elements with correspondent fluid rapidities 
Vi = Vi> see Fig. 8. In the case y\ ~ y\ the fluid elements are overlapping. The sizes I of these 
fluid elements are rather small, I ~ r« / ^* [SSI) as compared to the system's scales and also to the 





1 Except for an (artificial) situation with a box-like density profile 32 . 
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correspondent homogeneity lengths of direct pions. If the half-sum of pion transverse momenta is 
small, pt <C po, then the distance between the two fluid elements could be also small if y\ = y\ 
or fairly large, L 3> /, if y\ = —y\ ■ The former configuration corresponds to situations a) and 
b) in Fig. 8 at the moment of resonance decays r > tj . The later configuration c) increases 
significantly at small pt the long- interferometry radius, that corresponds to a Gaussian fit to the 
correlation function. It is enhanced additionally since in this case the velocities of the resonances 
are back-to-back directed (total momentum of pair is selected to be near zero), and they propagate 
for some time until they decay. At large pt > Po the distance between the two resonances is reduced 
according to the kinematics, especially when transverse flow takes place, which is why the resonance 
contribution to the interferometry volume dies out with pt much faster than their contribution to 
the spectra, see Fig. 6. 

Let us illustrate also the reduction of the ratio R ou t / Rside due to heavy resonance decays in the 
presence of transverse flows. For simplicity we omit the longitudinal coordinate and will keep in 
mind the linear with transverse radial flow and the Gaussian-like particle density profile at freeze- 
out. Let us consider at pt = Po the interferometry in transverse plane for the pions coming, as we 
described above, from the one sort of resonances characterized by m* ^> T and with pion decay 
momentum pq. In Fig. 9 we mark the four points a, b, c, d of emission of the resonances at freeze- 
out which are minimal and maximal in x (out-) and y (side-) directions at the given restriction 
for pion momenta. Those positions are dictated by the kinematics of the decay (vq = po/E n ) to 
guarantee the pion transverse momentum to be pt = Po and directed along the x axis. As one can 
see from Fig. 9, the maximal sideways extension R^ (along the y axis) of the emission region for 
the resonances producing the pions is equal to the maximal outward extension R ac of the region, 
and the maximal relative resonance velocities in the ^-direction Av x = 2vq is equal to the maximal 
relative velocities in the y-direction Av y = 2vq. If the density were constant, the effective emission 
sizes would be R x = R y . But since the density decreases with the radius, the effective outward 
size of the emission region, accounting for the probabilities of emission, will be less, R x < R y , since 
point R c is maximally remote from the centre as one can see from Fig. 9. Then we have also 
Av x < Av y effectively. It reduces the ratio of the interferometry out- to side- radii especially when 
mean life-time r* of the resonances are fairly large: R ou t = Rx + t*Av x < R s id e = Ry + t* Av y . 
The results of Fig. 7 demonstrate the total effect from the contributions of many resonance species 
in the 3D expanding hydrodynamic model. 

Note that the above effects are typical in the hydrodynamic picture with sudden freeze-out when 
the fraction of pions produced by resonances is fairly big as takes place in the case the chemical 
freeze-out. In the opposite case, mainly direct-direct and direct-resonance pion pairs are important 
for interferometry radii formation; the correspondent results are well known |29| I30| . 

3.3 Resonance decays and phase-space densities 

This detailed consideration allows us to clarify what we should measure to estimate the overpopu- 
lation of phase-space in the system at the latest stage of its evolution. First, the spectra of pions 
from decays of some fixed resonance species are primarily non-thermal; second, the resonances 
that produce the pion pair are emitted, in general, from different small elements of fluid separated 
by some distance that could be rather large at small transverse momenta of the pair. The pions 
of the pair are produced also at different (proper) times. This is completely unlike the emission 
of direct pions at freeze-out in the hydrodynamic picture: the pions with similar momenta p are 
emitted from the same fluid element, the value of momenta p dictates the size of that element 
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(the homogeneity length) and the HBT interferometry measures the correspondent size. Hence, 
the maximal phase-space density at each p is just the spectra n(p) divided by the interferometry 
volume and the APSD is less than it by a factor (2) 3 / 2 if one supposes a Gaussian phase-space 
density profile (at fixed p !) outside of that fluid element. Such a quantity, the APSD, for the case 
when x — p correlations are absent and the Gaussian profile is universal for all p was analysed in the 
pioneer Bertsch paper 1 . For pions created by the resonances there is no such transparent physical 
meaning of the ratio of the spectrum to interferometry volume, even if the resonance width is set 
to infinity and so they are only allowed to produce pions at the freeze-out hypersurface. Indeed, if 
we suppose the absence of transverse flows, then at small px of pairs there are, as was discussed 
above in detail (see also Fig. 8), the two rather small, I ~ Ty/T/m*, and well separated by distance 
L ~ 2rsinhr/i 3> I, fluid elements which emit the pions by means of resonances. Obviously, the 
maximal APSD here is proportional to n{p)/l and is a high value, while it will be estimated as 
a small value, ~ n(p)/L, according to the Bertsch prescription, when the interferometry radius, 
extracted by a Gaussian fit to the interferometry peak, is proportional to L . Of course, the 
correspondent correlation function in a large interval of qi is rather different than the Gaussian 
form; it is proportional to 1 + exp(— q 2 l 2 ) cos 2 (qL/2) , and this is the solution of the paradox - the 
simplest form, the Bertsch prescription, just fails here, but how does one distinguish such effects of 
the non-Gaussian contribution, when experimenters are fitting the total correlation function from 
many similar contributions to a Gaussian? The finite life-time of the resonances makes the situation 
even more complicated. We can conclude that the effective region (homogeneity volume), where 
the single pion spectrum from the resonance decays forms, is not associated with an interferometry 
volume of pion "resonance emission" that is topologically non-trivial and corresponds, in general, 
to the superposition of space-time separated Gaussian sources. In the next Section we propose a 
way to extract the actual APSD of pions in such a situation. 

4 Extraction of pion fugacities in A+A collisions. 

4.1 Sudden freeze-out vs continuous emission: how to treat the APSD? 

Our aim is to propose a simple method to estimate the overpopulation of phase-space that does 
not depend on the details of hydrodynamic evolution. Before discussing it, let us make some com- 
ments on the applicability of the hydrodynamic approach to the post hadronization stage in A+A 
collisions. As we discussed in Sec. 2, there is approximate equivalence between the chemically 
frozen hydrodynamics of a hadron gas and its evolution within the cascade approach in the tem- 
perature region of applicability of hydrodynamics. Out of this region, below the temperature 0.12 
GeV [S], the formation of spectra is basically continuous in time with fairly long "tails". That 
continuous emission is characterized by the emission function g(t, r, p) which, normally, is very 
far from thermal. The example of it is presented in Ref. [22] an d is based on the exact analytic 
solution Q of the Boltzmann equation for an expanding fireball. Unlike the very complicated 
structure of the emission function g(t, r, p) that was found in |22] . the interferometry radii, spectra 
and the APSD have simple and clear analytical forms and correspond to the thermal distribution 
function in the fireball before it starts to decay. It looks like there is some kind of duality [31] in 
the description of the spectra and the interferometry data based either on the (thermal) distribu- 
tion functions f(t,r, p) which characterize the system just before decay begins or on the emission 
function g(t, r, p) that describes the process of continuous particle emission after the system loses 
its simple thermal properties. The replacement of the complicated emission process by the simple 
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Landau j^B] prescription of sudden freeze-out at the lower boundary of the region of applicability 
of hydrodynamics, at T ~ 0.1 — 0.14 GeV (depending on the system size) is, of course, a rather 
rough approximation. Nevertheless, as we discussed above, the momentum-energy conservation 
laws, particle and entropy (partial) conservation as well as some symmetry properties minimize 
the correspondent deviation for observables. As for the average phase-space densities, which, as we 
discussed in Sec. 2, are approximately frozen during the hydrodynamic evolution, there is no reason 
to expect that they are changed significantly during the transition to their complete freeze-out (the 
APSD does not change during the stage of free streaming). 

Therefore, one can consider the analysis of the APSD in two aspects. The first one associates 
the observed APSD with the corresponding value in the thermal hadron medium just before its 
decay, that allows one to obtain, based on the distribution function f(t, r, p), the particle number 
and other densities in the medium at the pre-decay stage. On the other hand one can link the APSD 
with the properties of the particle continuous emission, and thus analyse the process of particle 
liberation described by the emission function g(t,r, p). If the decay of the medium starts at time 
tj, then one can compare the particle density just before decay n(tf, r) = f d 3 pf(tf,r, p) with the 
density of "emission capacity" of the source at point r, n em it(r) = J. dt J d s p g(t,r, p). Since 
the system continues to expand during the process of particle liberation, then n(tf,r) > n em it{r) 
and the equality is reached only in the case of a real sharp freeze-out when, in the non-covariant 
form, g(t, r, p) = f(t,r,p)5(t — tf). Thus the results of our analysis of the APSD based on the 
hydrodynamic picture with the Landau freeze-out prescription gives the possibility, on the one 
hand, to find the particle density in the thermal hadron medium just before its decay and, on 
the other hand, to make the upper estimate of the "emission capacity" density for the continuous 
process of particle liberation. 

It is clear, however, that the resonance decays at the post freeze-out stage destroy the above 
conception of duality in the APSD description since they seriously affect the post-hydrodynamic 
(pion) emission function g(t, r, p). The attempt to attribute the decays of resonances just to the 
sharp freeze-out hypersurface will make an estimate of the pion density to be artificially high. 
Indeed, the proper life-time 1/T of resonances averaged over all the species is nearly 2 fm/c that 
gives us, taking into account the Lorentz-dilation of life-times of resonances in the expanding swarm, 
an even larger value, say 2.5 fm/c, when the total number of resonances drops by a factor of 1/e. 
Since at the last quasi-inertial 3D-stage of expansion the system volume grows roughly proportional 
to t 3 , one can easily estimate, keeping in mind the typical (proper) time 10 fm/c for the freeze-out 
in A+A collisions at SPS and RHIC, that despite the resonance decays bringing about half of the 
pion yields the pion densities nevertheless decrease at the post freeze-out stage! It means that 
at the final decoupled stage of the system evolution the maximal densities, including the APSD 
of pions, are reached at the thermal freeze-out, at the boundary of the applicability region of the 
hydrodynamic approach. Thereby, taking into account also that resonance contributions to the 
APSD is physically obscure and their including complicates the treatment of the APSD very much 
(see the discussion at the end of Sec. 3) we will exclude the resonance contributions to the spectra 
and interferometry radii when estimating the maximal pion APSD. Since the effective temperatures 
of the total pion spectra and direct pions are approximately the same until rather high px , one 
can just reduce the spectra by the corresponding factor, 1/3 — 1/2, in accordance with the chemical 
freeze-out concept (HI IE]- As for the interferometry volume, the correction factor can be taken as 
0.6 for SPS energies and 0.5 for RHIC energies as follows from Fig. 4 and will be explained below. 
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4.2 A model independent method extracting fugacity 

Coming back to a hydrodynamic analysis of the spectra, interferometry radii and the APSD of 
thermal particles in Sec. 2, one can see (e.g. from Eqs. ()1U|) ) the simple link between the 
above values in the non-relativistic hydrodynamic models 

f ]3 AT 

U{t,p)) = -^/{^ 2 V mt ) (21) 

Eq. ()21j) expresses the basic Bertsch idea and can be generalized in covariant way for a relativistic 
situation |35| . If a system decays on some freeze-out hypersurface a, the PSD averaged over this 
hypersurface can be expressed through the two-particle correlation function C(p,q) of identical 
bosons as the following 

{fM) = Jf^(t,r,p)p^ = Wfl " q (22) 

Thus, to estimate the APSD one should use the experimental data on spectra and correlation 
functions or just the interferometry radii if C(p, q) — 1 is represented in standard Gaussian form, see 
H21j) . As for the left hand side of Eq. (I22|) , it is strongly model dependent since it is determined by 
unknown functions in space-time such as the fugacity £ = exp(/3/i) , temperature T = 1/(3 and flow 4- 
velocities u M in the local equilibrium Bose-Einstein distribution f(x,p) = (27r) _3 [£ _1 exp(f3p-u) — 1}. 
More concretely, it was demonstrated in Ref. ^3] that boost-invariant longitudinal flows in the 
system could bring significant error if the fugacity is extracted from (f(p)) assuming a static source 
2 . It was shown also that the extracted value of fugacity depends strongly on transverse flows - all 
of which makes estimates of this value essentially model dependent. Therefore the experimental 
analysis (see, e.g. j2Up of the overpopulation of phase space is carried out within some concrete 
model and, thus, depends on many assumptions such as distribution functions on a and the form 
of this freeze-out hypersurface in Minkowski space. 

Here we propose a method that allows one to extract the fugacity in a model independent 
way. Let us explain the basic idea using the oversimplified example when the freeze-out happens 
at uniform temperature and uniform particle densities across the whole freeze-out hypersurface 
enclosing the space-time of the expanding hydrodynamic system. Then the phase- space density 
averaged over the hypersurface and momentum, 

I^da,f(x,p) (2vr)3iV ' { } 

will be a suitable value which allows one to extract the fugacity at the freeze-out unambiguously, 
no matter to the form of the hypersurface and flow profile. Indeed, at each point x one can make 
the integration over p in the rest frame of the fluid element at the point x taking into account 
that in this comoving system the f(x,p) depends on the absolute value of p only (property of local 
equilibrium). Then, coming back to the common cm. system, and exploring the Lorentz- invariant 

2 Note, it is erroneously claimed in Ref. |14| that such a discrepancy vanishes at high pt, while the fugacity is 
still a factor of v2 larger (in the Boltzmann approximation) at high pr than one could estimate supposing that the 
source is static. 
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properties of the local equilibrium distribution function, one gets 



(2,f </(cr)> = /^MW ra^Hg (24 ) 
If, as we supposed, the temperature and particle densities are constant at freeze-out, then 



Sp 

(^exp^po)-!) 2 



(2tt)<« (/(a)) = 7 P , (25) 

J C- 1 cxp(/3p )-l 

and (/(cr)) = (/) e „ where (/) e „ is just the APSD for emission of a static thermal homogeneous 
source. 

Despite the extreme simplicity of this result, it cannot be used practically for two reasons. The 
first is that the integration in Eq. (J24*|) must be done over the whole closed freeze-out hypersurface, 
including, of course, butt-ends of the hydrodynamic tube. It implies the knowledge of one- and 
two-particle spectra in the whole momentum region (y,pr) which cannot be obtained by modern 
detectors. The second important reason why the formula (/(cr)) = (/) cannot be used is that 
the baryochemical potential and, hence, the temperature are not uniform across the longitudinal 
direction as the data on proton to antiproton ratios in different rapidity regions demonstrates [3^] ; 
see also the discussion about inhomogeneities of thermodynamic parameters in rapidity in Refs. 

0. Em 

Fortunately, a result similar to (|24[) can be obtained for the central rapidity interval where the 
spectra are usually measured and the thermodynamic parameters are nearly uniform. Let us start 
from the well known observation j^f)] that at the latest stage of a hydrodynamic evolution the 
longitudinal velocity distribution corresponds to the asymptotic quasi-inertial regime, vl = v x /t. 
It takes place even in the Landau model of complete stopping, and, of course, in the Bjorken 
model where this quasi-inertial regime exists at all times of hydrodynamic evolution. Further, 
if deviations of thermodynamic values such as temperature, baryochemical potential, etc. in the 
longitudinal direction near central fluid rapidity rj = are negligible within the maximal rapidity 

length of homogeneity [SHI; \J ~^T~ — ^> then au physical values at mid-rapidity can be calculated 
supposing boost-invariance. The latter implies that the distribution function which is applied for 
evaluation of the spectra at mid-rapidity near y = can be considered to depend only on the 
difference between fluid and particle rapidities, rj — y. As a result, the APSD at fixed particle 
rapidity (f(a,y)} is equal to the APSD (/(cr, 77)), 

f f 2 (t r nW^^E 

(/foiy)) = J ( / V l° , (26) 

//(t,r,p)^f 

at the fixed ratio r z /t = tanhry, corresponding to the fluid rapidity rj = y. Then using the Lorentz 
transformations as described above one can absorb all dynamic and kinematic characteristics into 
common factor V e ff(rj) = J ^| i w' t ( r i 7 ?) f° r numerator and denominator in Eq. (|2fij). Note, if the 
freeze-out hypersurface r(r) has non-space-like sectors, the measure of integration in (|26|) should 
be modified to exclude negative contributions to particle numbers at some momenta. As a result, 
it contains 9— functions like 9(p^ l n ll (x)) where is a 4-vector orthogonal to the hypersurface a 
[HHl E01- If a fluid element that crosses the time-like sector of the freeze-out hypersurface decays 
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preserving its total particle number, then the measure of integration is modified according to the 
prescription in Ref. [22]. In that case, after integration over all momenta at each point x the 
factorization of the numerator and denominator in (j2fi|) takes place again resulting finally in the 
common factor V e tf{rf) which cancels. Thereby, the expression (|26|) for the APSD can be used in 
the same form even when freeze-out hypersurface contains non-space-like sectors, if one considers 
the prescription from Ref. jHlj to describe the decay of the system on the time-like parts of the 
hypersurface. Finally, the common effective volumes are canceled and the result is 

f dp 

(2vr) U{a,y)) y=0 = dp (27) 

J C" 1 exp03po)-l 

Here the inverse of temperature (5 and fugacity £ characterize thermal properties of the decaying 
system at mid-rapidity and are supposed to be approximately constant within one unit of rapidity 
near rj = 0. On the other hand, the above value of the APSD is expressed through the transverse 
spectra and the interferometry radii at mid-rapidity similar to Eq. (|23|) . 



2-7T 5 / 2 J (^ Ro r sRl {^2-Km T dm T dy) 1 dnlT 
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Here we neglect interferometry cross-terms since they are usually rather small in the mid-rapidity 
region. The factor k, which will be calculated later from the results of the preceding Section, is 
introduced to eliminate contribution of short-lived resonances to the spectra and interferometry 
radii. It absorbs also the effect of the suppression of the correlation function due to long-lived 
resonances. 

Thereby we proved, using very natural assumptions, the possibility of extracting in a model in- 
dependent way the fugacity in expanding thermal systems at the stage of their freeze-out, realizing, 
thus, the program which was declared by Bertsch in his pioneer paper pQ. 

4.3 Pion APSDs and chemical potentials at SPS and RHIC 

Let us use our approach to find the chemical potential of direct pions at thermal freeze-out at SPS 
and RHIC. First we have to obtain the value of k in Eq. (|28|). This coefficient includes, in particular, 
the factor y/A (see, e.g., review j3J) arising because decays of the long-lived (I) resonances such as 
rj-, rj' -mesons give the contribution to the total pion spectra but not to the observed 

interferometry radii. The reason is that the (theoretical) correlation function due to pions produced 
by the long-lived resonances contains very narrow (much smaller than momentum resolution <7min 

of 

a detector) peak at q = on the top of smooth curve associated with bulk pion source. This narrow 
peak gives, in fact, no contribution to integral over q of (C(p,q) — 1) (see Eq. (|23[H in Bertsch's 
method, and, therefore, results just in reduction of the total value of the APSD. In accordance with 
the estimates made in Refs. |291 142) we choose the effective suppression factor A which measures 
the fraction of pion pairs containing no pions from long-lived sources, 



to be equal to 0.8; thereby k contains a factor yA ~ 0.9. Note, that experimental value of the 
suppression factor is, normally, smaller than the above theoretical value since the many effects, 
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except the decays of long-lived particles, could suppress the measured correlation functions. The 
most important among them are the single- and two-track resolution and particle misidentification. 
These effects are irrelevant to the physics of A+A collisions and we use, therefore, theoretical value 
of the suppression factor associated with long-lived resonances. 

Since the slopes of the complete pion spectra and thermal pions are very similar in a wide my 
region (see Fig. 4), we just reduce the experimental spectra by the ratio of the total pion number 
to the number of direct pions that follows from an analysis of the particle number ratios within 
the concept of chemical freeze-out. We use the reduction factor of 3 as is found in Ref. [S], also we 
estimate for comparison the fugacity supposing the correspondent reduction factor to be 2. The 
ratio of the complete to "direct" interferometry volumes can be fitted by a const/ ^Jm? function in 




the region shown in Fig. 4. We take into account this momentum dependence to define the effective 
correction factor to observed interferometry volume which we will accumulate again in k in Eq. 
(|28|) . The correspondent contribution is 1.65 for central Pb+Pb 158 AGeV collisions at CERN SPS 
and is 1.92 for Au+Au ^/sj^N = 130 GeV at BNL RHIC. Thus, k = 0.9 • (1/3) • (1.65 - 1.92) = 0.5 
(SPS) - 0.6 (RHIC). If one supposes the ratio of the total pion number to the direct one to be 2, it 
changes the values of k to k = 0.65 (SPS) - 0.7 (RHIC). We use also the latter values to estimate 
the sensitivity of the results to this ratio. 

To evaluate the APSD of direct pions (|28j) we utilize the following parametrization of the tt~ 
transverse spectra and interferometry radii. 



Rs = Cs,i/y/l + Cs^T and correspondent numerical parameters are taken from Ref. [33] (the 
CERES Collaboration). We use the approximation Ro = Rs until the minimal measured pt mo- 
mentum, pt = 0.125 GeV, and our analytical approximation of the CERES outward interferometry 
radii data for px > 0.125 GeV. 

For RHIC Au+Au = 130 GeV: 

The transverse spectrum is 2 nmfdL T dy = A ( ex P(- m T/T e f f ) - l)" 1 , T eff ss 0.218 GeV, and 
£~ ~ 249 (the STAR Collaboration, [IS]). We use here the Bose-Einstein parameterization of 
the transverse spectra from Ref. [3^1 since the integrated rapidity density of negative pions with 
this fitting function is closer to value £~ ~ 270 presented recently by the PHENIX Collaboration 
at midrapidity |46j than what one can get from the exponential parameterization of the transverse 
spectra, — g^- ~ 229 j3E|. The phenomenological parameterization of interferometry radii, Ri = 

Cl/^ttit, Rs = Cs,i/ \/l + Cs,2 + Cs^rnr, and the correspondent numerical parameters are taken 
from Ref. [37J of the PHENIX Collaboration. As one can see from Fig. 3 of Ref. [17], there is some 
discrepancy in the data on Rs radii between the STAR [3E] and PHENIX Collaborations (37J- To 
optimize uncertainties in the forthcoming estimates , we utilize for the parameter Cs,i the value 
8.75 fm that is the average between "PHENIX data motivated" value, 8.1 fm, and "STAR data 
motivated" value, 9.4 fm, [37J. We use the approximation Ro = Rs- 

The results of our calculations of the APSD according to 1)28(1 are collected in Table 2. The 
values found are used then to extract the pion fugacities based on our result (|27|). In Table 3 we 
present the pion chemical potentials and densities at different typical freeze-out temperatures. 

TABLE 2. The average phase-space densities (APSD) of negative pions for all particles ("raw") 

and for thermal ones. 
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without 
corrections 


K 

corrections 
for resonances 


(^r(f(a,y)) 
for thermal 
(direct) pions 


SPS 

^/s NN = 17.3 GeV 


0.210 


0.5-0.65 


0.105-0.137 


RHIC 

\J &nn = 130 GeV 


0.240 


0.6-0.7 


0.144 - 0.168 



TABLE 3. The chemical potentials, "raw" and for thermal pions, and thermal densities of 
negative pions vs critical densities extracted at two typical temperatures of the kinetic freeze-out. 





Tth MeV 


MeV 


MeV 


n n - fm~ 3 




freeze-out 


without 


for thermal 


thermal 




temperature 


corrections 


(direct) pions 


pion density 


SPS 


100 


78 


32-51 


0.014 - 0.017 


y/sNN = 17.3 GeV 








[0.065] 


SPS 


120 


70 


15 - 38 


0.022 - 0.027 


^/s NN = 17.3 GeV 








[0.095] 


RHIC 


100 


85 


54-64 


0.018 - 0.020 


\j &nn = 130 GeV 








[0.065] 


RHIC 


120 


78 


42-54 


0.028 - 0.031 


\l &nn = 130 GeV 








[0.095] 



All the values are related to the negative pions. The values in square brackets in the last 
column of Table 3 are the critical pion densities (in the rest frames of the fluid elements) at 
the correspondent temperatures for one (negative) component of the isospin pion triplet. As we 
can conclude, no exotic phenomena associated with the overpopulation of phase-space could be 
expected at the current energies. Some increase of the APSD takes place at RHIC; the pion 
chemical potentials are higher than at SPS and correspondingly the pion number density at RHIC 
is closer to the critical value for Bose-Einstein condensation. However the density values are still 
too small to expect the critical phenomena, such as serious spreading of multiplicity distributions, 
significant reduction of the interferometry radii, etc., which are associated with the Bose-Einstein 
condensation of pions in momentum space. Note, as was discussed above, our estimates of the pion 
densities at the end of the hydrodynamic stage of the evolution are, at the same time, the upper 
limits for the "emission capacity" density of the source of the direct pions. 

The APSDs per pion component in the chemically equilibrated neutral pion gas are 0.07 and 
0.09 for T th = 0.100 GeV and T th = 0.120 GeV respectively. As one can see from the Table 2, the 
estimated APSD for thermal pions are somewhat higher than in the case of chemical equilibrium. 
Thereby some corrections to observables conditioned by the partial overpopulation of the pion 
phase-space should be taken into account in advanced models of A+A collisions. 

5 Conclusions 

An approach which allows one to estimate the chemical potential and, thus, the overpopulation 
of particle phase-space at thermal freeze-out in a way that does not depend on the flow profile 
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and form of the freeze-out hypersurface has been developed. It realizes Bertsch's basic idea to 
find from the interferometry data the phase-space density averaged over both configuration and 
momentum spaces and compare it with the correspondent value for static source with equilibrium 
(homogeneous) Bose-Einstein thermal distribution to estimate thereby the overpopulation of the 
phase-space. We have proved that under the assumption of uniformity of particle density at the 
freeze-out hypersurface the flows and form of the hypersurface are absorbed into a factor that is 
canceled when the APSD is determined. It is worth to note that such a "cancellation" takes place 
not only for the APSD. For example, as one can find easily under the same conditions, the thermal 
pion entropy per particle also does not depend on the form of the hypersurface and flows. 

Another problem that has been analysed is the resonance contributions to pion phase space 
densities. Because of chemical freeze-out an essential fraction of pions, about half, will be produced 
by short-lived resonances after thermal freeze-out. It is found that the decays of many heavy 
resonances correct significantly the interferometry radii of direct pions. It could result in up to 
a 50 per cent increase of the correspondent interferometry volume at small transverse momenta. 
Since the heavy resonances have, typically, large momentum po of pions produced, the resonances 
contribute to the interferometry radii up to pr ~ 1 GeV. One of the important consequences of this 
is the ratio of the outward interferometry radius to the sideward one, Rout/ Rside, that becomes less 
than unity due to heavy resonance decays in contrast to standard results R ou t / R s ide > 1 for direct 
pions in the hydrodynamic models (except for the so-called "blast- wave" model [31], see Footnote 
1). It could be a step towards understanding the HBT puzzle. 

The correlation function of two identical pions arising from decays of heavy resonances of 
fixed species, even if one neglects their life-times, is essentially non-Gaussian, especially at small 
pion transverse momenta. Then, if one gets the interferometry radii from a Gaussian fit to the 
correlation function in the region of the interferometry peak, the longitudinal radii will be larger 
than the homogeneity lengths of regions where the spectra of "resonance pions" are formed. This 
violates the natural physical correlation (f(t,p)) ~ ^p^/Vint that is typical for direct particles. 
Accounting for the finite life-times of resonances makes the picture even more complicated: the 
emission is now spread out in time direction and therefore the Vi n t cannot be attributed to the 
thermal freeze-out hypersurface. We have argued thereby that resonance contributions should be 
excluded when the averaged pion phase-space densities at the freeze-out stage are determined. Our 
estimates of the total and thermal spectra and the ratios of correspondent interferometry volumes at 
different intensities of transverse flow, including what are typical for SPS and RHIC energies, gives 
the possibility of restoring the APSD of thermal pions at thermal freeze-out. We argue that the 
values found are the maximal (total) pion APSD at and after the decoupling of the hydrodynamic 
system. 

We found that the overpopulation of thermal pions at the thermal freeze-out, 100 < Tth < 120 
MeV, is associated with the pion chemical potentials 15 — 51 MeV for Pb+Pb(Au) 158 AGeV 
collisions at CERN SPS and 42 - 64 MeV for Au+Au ^v]v = 130 GeV collisions at BNL RHIC. 
The estimates of the correspondent values in the recent paper |49| based on the concrete hydro- 
dynamic parametrization are close to our results for SPS energies, while for RHIC energies the 
pions chemical potentials estimated in Ref. [IS] are somewhat higher. The thermal pion densities 
turn out to be essentially smaller than the critical ones. The correspondent values at the SPS 
energies are 0.014 - 0.017 < 0.065 fm~ 3 for T th = 100 MeV and 0.022 - 0.027 < 0.095 fur 3 for 
T th = 120 MeV. For the RHIC energies we found 0.018-0.020 < 0.065 fm~ 3 for T th = 100 MeV and 
0.028 - 0.031 < 0.095 fm~ 3 for T th = 120 MeV. Thereby the multibosonic effects at those energies 
would be considered rather as a correction factor than as an important physical phenomenon. 
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It looks puzzling that the averaged phase- space densities of thermal pions at RHIC are higher 
than at SPS, 0.144 — 0.168 > 0.105 — 0.137 while the particle freeze-out densities are comparable. 
The detailed analysis of the evolution of the APSD could explain such a behaviour. Indeed, if 
one imagines that at some initial time hadronic system has no significant transverse flow, then the 
averaged PSD is proportional to the PSD and particle density, as in a hot, dense gas in a box, 
and all the densities are high. However, when the system expands and transverse flow develops, 
the behaviors of the APSD and particle density are completely different. While the pion density 
is reduced dramatically, at the last stage as i -3 , the pion APSD changes (it can even increase) 
rather slowly. Thereby the systems reach kinetic freeze-out at some typical particle densities that 
are similar for different collision energies, and, at the same time, their APSD could increase with 
energy since it "memorizes" the higher initial hadronic densities at higher energies. It is found 
that such a "memory" is conditioned by the conservation of entropy and particle numbers during 
chemically frozen hydrodynamic evolution of the mixture of hadronic gases. 

The main aspect of the HBT puzzle concerns the energy dependence of the interferometry 
volume that turns out essentially smaller than expected at the RHIC energies. Our results could 
shed light on this phenomenon. Since the pion APSD and effective transverse temperature do not 
change essentially during the evolution, the "interferometry volume", Vi n t ~ (dN j dy) / {T^j t (/(£))), 
defined formally at any time during the evolution of hadronic system, also changes rather slowly 
with time. The concrete mechanism of the effect is that, if the intensity of the transverse flow grows 
during the evolution, the contribution of an enlarged geometric radius of the expanding system to 
the interferometry volume is almost compensated by a stronger reduction of the interferometry 
radii due to the strengthening of the flow (see, e.g., [23]). The initial effective volume is, thus, 
partially frozen in the form of the interferometry volume. There is a rough similarity between the 
effective, or "interferometry" volumes at higher (RHIC) and smaller (SPS) energies because the 
larger configuration volume at hadronization stage at RHIC (particle densities at the hadronization 
are approximately the same as at the SPS) is partially compensated by the higher intensities of 
flow. Therefore, the unexpectedly small increase of the interferometry volume at RHIC is caused 
by strengthening of transverse flows at higher energies. During the evolution of hadronic system 
the APSD is reduced rather slowly ("frozen") while the particle density falls down quickly until 
it reaches the standard conditions for system decay. That is the reason why the observed value 
of the APSD has no direct link to the freeze-out criteria and final thermodynamic parameters, 
being connected rather to the initial phase-space density of hadronic matter formed in relativistic 
nucleus-nucleus collisions. 
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Figure 1: Evolution of the pion effective temperature T e ff (solid line) versus the temperature of the 
system T (dashed line) in the non-relativistic one-component Bj or ken-like model with transverse 
expansion. The plots correspond to initial conditions: To = 0.1 GeV, time to = 5 fm/c, transverse 

radius Rq = 5 fm and initial transverse flow R (to) = 0.3. 
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Figure 2: The average phase-space densities {f(r,p)) at typical proper times: at hadronization, 
r = 7.24 fm/c, (solid line) and at kinetic freeze-out r = 8.9 fm/c (dashed line). The dot-dashed line 
corresponds to the "asymptotic" time r = 15 fm/c of hydrodynamic solutions continued beyond 
kinetic freeze out. The initial conditions of hydrodynamic expansion are taken from Table 1 and 
based on Ref. [S]. 
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Figure 3: The space densities n(r) taken at typical proper times: at hadronization, r = 7.24 fm/c, 
(solid line) and at kinetic freeze-out r = 8.9 fm/c (dashed line). The dot-dashed line corresponds 
to the "asymptotic" time r = 15 fm/c of hydrodynamic evolution continued beyond kinetic freeze 
out. The initial conditions of the hydrodynamic expansion are taken from Table 1 and based on 
Ref. 0. 
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Figure 4: The pion transverse spectra for the complete pion yields (solid lines) and for "direct" pions 
only (dashed line) in the cases with and without transverse flow: vr = 0, 0.4 and 0.6. The shape 
of the pion spectra from the fixed heavy resonance species /2(1270) is presented for illustration by 
dot-dashed line for vr = 0. The overall normalization is arbitrary. 
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Figure 5: The pion transverse spectra for the complete pion yields (solid lines) and for "direct" 
pions only (dashed line) taken at a few typical temperatures without transverse flow, vr = 0. The 
overall normalization is arbitrary. 
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Figure 6: The ratio of the complete interferometry volume Vtot, that includes effects of the res- 
onance decays to pions, to the interferometry volume for "direct" pions taken at a few typical 
intensities of flow. The resonance contributions are calculated in accordance with the chemical 
freeze-out conception, when approximately 2/3 of pions are produced by resonance decays. The 
dashed line corresponds to vr = 0, solid line: vr = 0.4 and dash-dotted line: vr = 0.6. The other 
parameters are taken from Ref. [S] as described in the text. 
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Figure 7: The R ou t to R s id e ratio for complete interferometry radii (solid line) vs ones for "direct" 
pions (dashed line) at transverse flow vr = 0.6 which is associated with RHIC energies. 
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Figure 8: The two- pion spectra formation at pr — and y ~ due to decays of heavy resonances 
involved in the longitudinal flow before and at the freeze-out time t/. g ., see the text for details. 
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Figure 9: The formation of the pion transverse spectra for decays of the heavy resonances involved 
in transverse flows as it is described in the text. The solid arrows are associated with radial flow, 
the dashed arrows correspond to emission of pions from the resonances. 
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